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Abstract 

We update our approximate parametrizations of the three-loop sphtting functions for 
the evolution of unpolarized parton densities in perturbative QCD. The new information 
taken into account is given by the additional Mellin moments recently calculated by 
Retey and Vermaseren. The inclusion of these constraints reduces the uncertainties of 
our approximations considerably and extends their region of applicability by about one 
order of magnitude to lower momentum fractions x. 



PACS: 12.38.Bx, 13.60.Hb 



In order to achieve a high accuracy of the predictions of perturbative QCD for hard 
processes, the calculations need to transcend the standard next-to-leading order (NLO) 
approximation. For processes with initial-state hadrons, the next-to-next-to-leading order 
(NNLO) expressions include the three-loop splitting functions. The computation of these 
functions is under way but will not be completed in the near future 0. 

Partial results have already been obtained [3-9], however, most notably the five lowest 
even-integer moments for the flavour non-singlet combination entering electromagnetic 
deep-inelastic scattering (DIS) 0, and four moments for the singlet splitting functions 
In refs. JTO], |ll|] we have demonstrated that this information — due to the smoothening 
effect of the ubiquitous convolution with the initial parton densities and the small size 
of the corrections — is fully sufficient for momentum fractions a; > 0.1 and leaves only 
small uncertainties down to x ~ 10"'^ at scales above about 10 GeV^. We have provided 
approximate parametrizations for the three-loop splitting functions, including quantitative 



estimates of their residual uncertainties. These results have already been applied ||T2[ to 



structure functions in DIS and Drell-Yan cross sections at hadron colliders, for which the 



subprocess cross sections have been computed up to NNLO [O P 



Very recently the fixed-moment calculations of refs. 0, ^ have been extended, using 
improved computing resources, up to the twelfth moment 0. For the first time also (odd) 
moments of the three-loop valence splitting functions have been obtained there. These 
results provide a severe check of our approximation procedure, as the latter led to rather 
tight predictions for the (tenth and) twelfth moments. This test is passed by the results 
of refs. |1T0| , p!T[] . In this letter we update these parametrizations by including the moments 
of ref. 0] in the derivation. As a result the residual uncertainties are greatly reduced, 
and the region of safe applicability is extended by about one order of magnitude in x, an 
improvement most relevant for applications to structure functions at HERA [p^ . 

Our notations for the parton densities and splitting functions are as follows: the non- 
singlet combinations of quark and antiquark densities, and qi, are given by 

Nf 

?NS = Qi ± Qi - {qk ± Qk) , ^NS = J2(^r-qr) ■ (1) 

r=l 

Nf stands for the number of effectively massless flavours. The corresponding splitting 
functions are denoted by P^g and -Pns = -Pns + -^ns- The latter function, P^g, occurs for 
the first time at O(a^). Except for the vanishing of the first moment, it was unknown 
before the calculation of ref. . The evolution equations in the singlet sector are written 

d ' S 1 _ I -Pns + -Pps Pqg i ~ / S 



Here g represents the gluon density, and ® denotes the Mellin convolution. The expansion 
of all these splitting functions in powers of the running coupling constant reads 

Fix, as) = as P^^\x) + al P^^\x) + al P^^^\x) + . . . with a, = a,/47r , (3) 
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if the renormalization and factorization scales are identified, Hr = /^/- The additional 
terms for fir ^ fJ'f are exactly known up to NNLO and need not to be considered here. 

Let us briefly illustrate our approximation procedure for the case of the Nf term, 
Pgg\{x), of the gluon-quark splitting function Pgg which dominates the small- a; evolution 
of the quark densities. Taking into account the small-x result of ref. 0, the expected 
form of this function in the MS scheme employed throughout this paper is given by 

Pifii^) = jl^m ln"^(l -x) + UooM + j2Bn\n^x+--^— , (4) 
^1 ^1 X 27 X 

where /smooth is finite for < x < 1. We choose three or two of the large-x logarithms, 
one or two smooth functions (mainly low powers or simple polynomials of x) and two of 
the small-x terms {x~^ and Inx or In^ x). Their coefficients are then determined from the 
known six moments [§, Varying these choices we arrive at the about 50 approximations 
compared in Fig. 1. The two representatives spanning the error band for most of the x- 

(2) 

range are finally selected as our estimates for Pgg'i and its residual uncertainty. 

Analogous procedures are applied to the N'j and Nj terms of all other functions 
P^^\x). The non-singlet N'j contribution is known 0. The singlet Nj pieces are smaller 
in absolute size and uncertainty than the Nj and Nj terms, hence for them it suffices 
to select just one central representative. Note that the previous information, four (five) 
moments in the singlet (-Pns) sector, respectively, was not sufficient to fix the coefficients 
of the subleading small-x terms oc x~^ (In^ x) from the moments. Thus we had to resort 
to conservatively varied, educated guesses inspired by the NLO results. Except for the 
case of p(2)+ (where the rightmost pole in Mellin space is not one, but two units away 
from the lowest calculated moment) we can now dispense with these ad hoc estimates. 

We now write down our improved parametrizations, using the abbreviations 

Lo = Inx , Li = ln(l — x) . (5) 

Two approximations, denoted by P4 and -Pg , are provided for each function. Where 
both are present, the Nj and Nj terms have been combined such that the error bands 
are maximized at small x. The averages 1/2 [A + i?] represent the central results. 

Our new expressions for the non-singlet splitting functions P^^ read 

pg^(a;) = 1185.229 (1 - x);^ + 1365.458 5(1 - x) - 157.387 - 2741.42 x^ 

- 490.43 (1 - x) + 67.00 + 10.005 Ll + 1.432 

+ Nf {-184.765 (1 - x);^ - 184.289 5(1 - x) + 17.989 + 355.636 x^ 

- 73.407 (l-x)Li + 11.491 + 1.928 Ljj} + PisU^;) 
pg-(a;) = 1174.348 (l-x);^ + 1286.799 5(1 -x) + 115.099 + 1581.05 Ll 

+ 267.33 (1 - x) - 127.65 Ll - 25.22 + 1.432 
+ Nf {-183.718 (1 - x);^ - 177.762 5(1 - x) + 11.999 Ll + 397.546 x^ 

+ 41.949 (1 - x) - 1.477 L^- 0.538 + P}§2ix) (6) 
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and 

PjJ|+ (x) = 1183.762 (1 - x)-^ + 1347.032 5(1 - x) + 1047.590 Li - 843.884 x^ 

- 98.65 (1 - x) - 33.71 + 1.580 (L^ + All) 

+ iV/ {-183.148 (1 - x)-^ - 174.402 5(1 - x) + 9.649 + 406.171 x^ 
+ 32.218 (l-x) + 5.976 + 1.60 Lj^} + Piiya^) 
pg+ (x) = 1182.774 (1 - + 1351.088 5(1 - a;) - 147.692 Ll - 2602.738 

- 170.11 + 148.47 Lo + 1.580 (L^ -All) 

+ AT/ {-183.931 (1 - x)-^ - 178.208 5(1 - x) - 89.941 Li + 218.482 

+ 9.623 + 0.910 - 1.60 L|^} + Piiiyx) . (7) 
The In^x terms in Eqs. (P) and (^ stem from ref. 0]. The exactly known N'j contribution 

-Pns 2(^) both cases is given in Eq. (4.13) of ref. [|10|- Pns {x) and Pjljg (x) are 
compared at x < 1 in Fig. 2, for A^j = 4, to our previous approximations based on one 
moment less for Pj^g and mostly indirect information on Pj^jg (x). Our present results 
are consistent with, but considerably more accurate than those of ref. [|TU . 



In contrast to Pqg and P^g, the transition from one to two loops leads only to ln^(l — x) 
terms in Pns and Pgg. Assuming correspondingly that no large-x logarithms beyond 
In'fl — x) occur in Pj^g* , the coefficients of the 1/(1 — x)_|- term of the (also about 50) test 

^1 



functions considered for P^g* (x) cover the range 1167.3 . . . 1190.3 for the A^9 part, and 



— 184.8 . . . —183.1 for the Nj contribution. The findings for PNg*^, where one moment 
less is known, are consistent with these results. Note that terms [ln'^(l — x)/(l — x)] + 
with > 1 do not occur in the MS splitting functions, as proven in ref. [O. 



The difference Pjfg = P^g — P^g and the pure-singlet splitting function Ppg in Eq. (0) 
result as Nf{P^.^^—P^.^J and Nf{P^.g^+P^^gJ , respectively, from the flavour independent 
('sea') parts of the quark-quark and quark-antiquark splitting functions. The former 
combination carries the colour factor d"''"^d"'^'^ which does not occur in P^c ; this fact 



facilitates the separations of the two terms in the results of ref. [H. Both Pns(x) and 



NS ' 

Ppg(x) vanish at x = 1, but are large at small x. The parametrizations selected for the 
contributions to these two functions are given by 

^NS,1(^) = A^/{(l-a;)(-1441.57x2 + 12603.59 X- 15450.01) + 7876.93 xLg 

- 4260.29 Lo - 229.27 Ll + 4.4075 L^} 
^NsS(^) = A^/ {(1 -x)(-704.67x^ + 3310.32x2 + 2144.81 X- 244.68) 

+ 4490.81 x^Lo + 42.875 Lq - 11.0165 L^} (8) 

pg^(x) = A^/ {(l-x)(-229.497Li- 722.99x2 + 2678.77- 560.20 x"^) 

+ 2008.61 Lo + 998.15 Ll - 3584/27 x^^Lq} + P^g^a;) 
pg^(a;) = A^j {(1 -x)(73.845L2 + 305.988 Li + 2063.19 X- 387.95 x"^) 

+ 1999.35 xLo- 732.68 Lo- 3584/27 x^^Lo} + P^%ix) (9) 



and 
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with 

P^^l^{x) = ATj {(l-a;)(-7.282Li- 38.779x2 + 32.022 X- 6.252 + 1.767 s"^) 

+ 7.453 L^} . (10) 

The {lnx)/x term in Eq. (|^) has been derived in ref. 

Our new approximations for the off-diagonal singlet splitting functions Pj^^^ and Pj^'^^ 
are given by 

P^I]a{x) = Nf {-31.830 Ll + 1252.267 Li + 1999.89 x + 1722.47 + 1223.43 Ll 



1334.61 a;-^ - 896/3 x-^Lo} + P, 



(2) 
99,2 



X) 



P^feix) = Nf {19.428 + 159.833 Lf + 309.384 Lf + 2631.00 (1 - x 

3(2) 
99,2 



67.25 Ll - 776.793 x'^ - 896/3 x-^Lq} + P^%{x) (H) 



with 



Pql^x) = A^^ {-0.9085 L^- 35.803 Ll- 128.023 + 200.929 (1 -x) 

+ 40.542 Lo + 3.284 x"^} , (12) 



and 



P^f j^{x) = 13.1212 L^ + 126.665 L? + 308.536 L? + 361.21 - 2113.45 Lo 

- 17.965 x-^Lo + A^/ {2.4427 L^ + 27.763 L? + 80.548^2 

- 227.135 - 151.04 L2 + 65.91 x~^Lo} + PgJ.Ua;) 

pj2)^(x) = -4.5108 Lf- 66.618 L?- 231.535 L^- 1224.22 (l-x) + 240.08 
+ 379.60 x-\Lo + 4) + Ar/{-1.4028 - 11.638 + 164.963 Li 

- 1066.78 (l-x) - 182.08 L2 + 138.54 x-i(Lo + 2)} + P<S5(^) (^3) 

with 

Pgfd^) = A^J {1.9361 L? + 11.178 Li + 11.632 - 15.145 (l-x)+ 3.354 Lo 

- 2.133 x"^} . (14) 

Unlike the case of P^g^ discussed above, the coefficients of the leading small-x terms 
(lnx)/x have not been derived for Pj^^ up to now. Thus these coefficients in Eq. 
have been determined, as before ilTlIl, from the available moments. 



The 1/[1 — x]4. soft-gluon contributions to the one- and two-loop gluon-gluon splitting 
functions, Pj°^ and Pgl\ are related to their quark-quark (non-singlet) counterparts by a 
factor Ca/Cp- The same holds for the iV| terms at third order Assuming that this 

relation holds generally for Pgf ■, the results given below Eq. (|^) can be employed. In this 
way we arrive at the approximate expressions 
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P^l]^{x) = 2626.38 (1 - x);^ + 4424.168 5(1 - x) - 732.715 Lj - 20640.069 x 

- 15428.58 (1 - x^) - 15213.60 + 16700.88 x'^ + 2675.85 x^^Lq 
+ Nf {-415.71 (1 - x)-^ - 548.569 5(1 - x) - 425.708 Li + 914.548 x^ 

- 1122.86 -444.21^2 + 376.98x^1 + 157.18 x^^Lo} + Pgf,2ix) 
P^f j^{x) = 2678.22 (l-x);^+ 4590.570 5(1 -x)+ 3748.934 Li + 60879.62 X 

- 35974.45 (1 + x^) + 2002.96 + 9762.09 x"^ + 2675.85 x'^Lq 
+ Nf {-412.00 (1 - x);^ - 534.951 5(1 - x) + 62.630 + 801.90 

+ 1891.40 Lo + 813.78 Lg + 1.360 x"^ + 157.18 x"^Lo} + Pgf,2{x) (15) 

with 

PglU^) = iV/ {-16/9 (1 -x);^ + 6.4882 5(1 -x) +37.6417x2 - 72.926 X 

+ 32.349 - 0.991 ^2 + 2.818 x"^} . (16) 

The (lnx)/x terms in Eq. (P!3| ) have been determined in ref. in a scheme equivalent to 
the DIS scheme up to NNLO. The transformation to MS can be found in ref. ||11|| . 



The uncertainty bands for the three-loop singlet splitting functions resulting from 
Eqs. (^-(0) are displayed in Fig. 3 for Nf = 4. As in the non-singlet cases consid- 
ered above, our new parametrizations considerably improve on the previous uncertainties. 
While the small-x behaviour of our approximations obviously depends on the results of 
refs. H, it is worthwhile to note that reducing the Nf coefficient of (lnx)/x in Pj^) 
a factor of two does not lead to approximations outside the error band in Fig. 3. 

In Fig. 4 we finally illustrate the impact of the NNLO terms on the scale derivatives (0) 
of the singlet quark and gluon densities. As in ref. Illf] the initial conditions are chosen as 



xS(x,/iJ) = 0.6x"°-3(l -x)^-^(l + 5x°-^) 

xg{x,ij}) = 1.0x-°-3'(l-x)^ (17) 



and 



«,(/i2 = fx}) = 0.2 , (18) 

corresponding to fij- ~ 30 GeV^. Under these conditions the residual uncertainties of the 
three-loop contributions amount to about ±2% or less down to x ~ 10"'*, even if the bands 
in Fig. 4 were increased by 50% in order to account for any possible underestimate of the 
errors. At lower scales the fiatter small-x shapes of the quark and gluon densities, together 
with the larger a^, lead to larger uncertainties at small x. At x = 10~^ and /ij ~ 3 GeV^, 
for example, they reach about ±4% and ±3% for the singlet quark and gluon derivatives, 
respectively, for standard NLO distributions like CTEQ4M |jl6| . These numbers represent 



improvements by about a factor of three on our previous results |11]. Thus our present 



approximations, based on the results of ref. 0, facilitate a reliable NNLO evolution of 
unpolarized parton densities down to, at least, /ij > 10 GeV^ and x > 10""^. 

Fortran subroutines of the above approximations of the three-loop splitting functions 
can be obtained via email to neerven@lorentz.leidenuniv.nl or avogt@lorentz.leidenuniv.nl. 
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Figure 1: Approximations of the Nj part Pqg\ of tlie tliree-loop splitting function P^g\x), 
as obtained from the six lowest even-integer moments 0, ^ together with the leading 
small-x term of ref . . The full curves represent those functions selected for Eq. dlTD . 
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Figure 2: Top: Our new approximations of Pj^s (^) ~ 4, as obtained from 

Eq. togetlier with Eq. (4.13) of ref. |T0|. The dotted curves represent our previous 
parametrizations ||TD[. Bottom: The same for Pj^g (x) using Eq. (|^). 
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Figure 3: Our new approximations of the singlet splitting functions P'>'^\x) for Nf = 4 



Pjg^ is obtained by adding P^s~^{x) of Fig. 2 (separately shown by the dash-dotted curve) 



and Pf^ix) of Eqs. (D and (|T0|). Our previous parametrizations [O] are also shown. 
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Figure 4: The size and remaining uncertainties of the NNLO corrections for the scale 
derivatives, S = dT, / din fi'j and g = dg/dln^'j, of the singlet quark and gluon densities 
at = /i^ ~ 30 GeV^ (a^ = 0.2). The input densities are specified in Eq. (|T^. 
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